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ABSTRACT 

Several statistics are applied to groups and galaxies in groups in the Two degree Field 
Galaxy Redshift Survey. Firstly we estimate the luminosity functions for different 
subsets of galaxies in groups. The results are well fitted by a Schechter function with 
parameters M* - 51og(/i) = -19.90 ± 0.03 and a = -1.13 ± 0.02 for all galaxies in 
groups, which is quite consistent with the results by Norberg et al. for field galaxies. 
When considering the four different spectral types defined by Madgwick et al. we 
find that the characteristic magnitude is typically brighter than in the field. We also 
observe a steeper value, a = —0.76 ± 0.03, of the faint end slope for low star-forming 
galaxies when compared with the corresponding field value. This steepening is more 
conspicuous, a = —1.10 ± 0.06, for those galaxies in more massive groups (M. :> 
10 14 /i _1 M Q ) than the obtained in the lower mass subset, a = —0.71 ± 0.04 (Ai < 
IO^^Mq). 

Secondly, we compute group total luminosities using Moore, Frenk & White pre- 
scriptions. We define a flux-limited group sample using a new statistical tool devel- 
oped by Rauzy. The resulting group sample is used to determine the group luminosity 
function finding a good agreement with previous determinations and semianalytical 
models. 

Finally, the group mass function for the flux-limited sample is derived. An excellent 
agreement is obtained when comparing our determination with analytical predictions 
over two orders of magnitude in mass. 

Key words: galaxies: clusters: general - galaxies: luminosity function, mass function- 
galaxies: statistics. 



1 INTRODUCTION 

Groups of galaxies constitute one of the most suitable lab- 
oratories for the study of properties of intermediate galaxy 
density environments and their consequences on the pro- 
cess of galaxy formation and evolution. Furthermore, sev- 
eral hints about the large scale structure of the universe and 
how structures evolve in the universe can be drawn from the 
statistical studies of groups and their properties. Some of 
them, for instance the luminosities and morphological types 
of their member galaxies, are sensitive to the processes of 
mergers and interactions between individual galaxies inside 
a potential well. Meanwhile, other properties as group abun- 
dance as a function of total luminosity or mass, can provide 
constraints to hierarchical clustering scenarios and cosmo- 
logical models. In this paper, we provide a robust statistical 
study about the luminosities of group galaxy members, and 
also on global group properties such as total luminosities 



and masses, using one of the largest group catalogues at the 
present. The accuracy of these determinations allow us a fair 
comparison with analytical and semianalytical predictions. 

Most studies on the environmental dependence of the 
galaxy luminosity function have been carried out either in 
the field or in rich clusters using the largest two dimensional 
catalogues and small redshift surveys. Regarding field galaxy 
luminosity functions, several works have been devoted to 
its determination in the last decade (Loveday et al. 1992, 
Marzke, Huchra & Geller 1994, Lin et al. 1996, Zucca et 
al. 1997, Ratcliffe et al. 1998). With the advent of large 
redshift surveys such as the Two degree Field Galaxy Red- 
shift Survey (2dFGRS) and the Sloan Digital Sky Survey 
(SDSS), more reliable statistical results have been obtained. 
The field luminosity functions determined by Blanton et al. 
(2001) (SDSS) and Norberg et al. (2002) (2dFGRS) show 
an excellent agreement, finding a luminosity function accu- 
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rately described by a Schechter function with parameters 
M 6 * 7 - 51og(/i) ~ -19.66 and a ~ -1.21. In particular, the 
2dFGRS allowed the determination of the field luminosity 
function for galaxies of different spectral types (Madgwick 
et al 2002) finding a systematic steepening of the faint end 
slope moving from passive (a = —0.54) to active (a = —1.5) 
star forming galaxies, and also a corresponding faintening of 
M* . 

A controversial issue about the luminosity function of 
galaxies in clusters is the steepening at the faint end, com- 
pared with field galaxy luminosity function (Valotto et al 
1997, Trentham 1997, Lopez-Cruz et al 1997). Recently Goto 
et al. (2002) computed a composed luminosity function us- 
ing 204 clusters taken from SDSS (York et al. 2000). They 
found that the slopes of the LF's become flatter toward red- 
der color band and that have brighter characteristic magni- 
tude and flatter slopes than the field LF. 

The lack of statistical results on smaller overdensities 
such as groups of galaxies is mainly due to the fact that two 
dimensional identification privileges the largest overdensi- 
ties. Analyzing a sample of 66 groups of galaxies identified 
in redshift space, Muriel, Valotto & Lambas (1998) find a 
fiat faint end for the galaxy luminosity function in groups 
(a ~ —1.0) compared with the luminosity function in clus- 
ters where a large relative number of faint galaxies is present. 
It is important to remark that most of the luminosity func- 
tion estimations in groups and clusters of galaxies are com- 
puted subtracting background and foreground contamina- 
tion due to the lack of spectroscopic information for galaxy 
members. 

One step further, in order to understand the transition 
between galaxy and galaxy systems luminosities, is the com- 
putation of the luminosity function of galaxy groups. Moore, 
Frenk & White (1993), reanalysing the groups in the Center 
for Astrophysics (CfA) redshift survey, developed a method 
for the estimation of the total luminosity of groups identi- 
fied in magnitude-limited galaxy surveys. This method al- 
lowed them the computation of the luminosity function of 
galaxy systems for a sample of 163 groups with at least 
tree members. Another attempt to determine the luminos- 
ity function of virialized systems was made by Marinoni, 
Hudson & Giuricin (2002). They used the Nearby Optical 
Galaxy (NOG) sample, which comprise ~ 7000 galaxies with 
cz < 6000km s~ x and B < 14, finding a very good agree- 
ment with Moore, Frenk & White (1993) previous determi- 
nation. 

On the other hand, the abundance of haloes as a func- 
tion of mass constitutes a key point in both, the determina- 
tion of a cosmological model and the understanding of the 
structure collapse. At the present, the more popular mod- 
els for halo abundance are the analytical model of Press & 
Schechter (1974) for spherical collapse, the Sheth & Tor- 
men (1999) model for ellipsoidal collapse and Jenkins et al. 
(2001) fit obtained from numerical simulations. Many efforts 
have been made to determine the mass function of galaxy 
systems from observations (Bahcall & Cen 1993, Biviano et 
al. 1993, Girardi et al. 1998). Recently, Girardi & Giuricin 
(2000) have computed the mass function for a sample of 
nearby loose groups by Garcia (1993). They found the group 
mass function to be a smooth extrapolation of the cluster 
mass function and a reasonable agreement with the Press & 
Schechter (1974) predictions. 



Currently, one of the largest group catalogue was con- 
structed by Merchan & Zandivarez (2002). They have iden- 
tified groups on the 2dF public 100K data release using a 
modified Huchra & Geller (1982) group finding algorithm 
that takes into account 2dF magnitude limit and redshift 
completeness masks. This catalogue constitutes a large and 
suitable sample for both, the study of processes in group en- 
vironment and the properties of the group population itself. 
The global effects of group environment on star formation 
was analysed by Martinez et al. (2002) using this catalogue. 
They have found a strong correlation between the relative 
fraction of different galaxy types and the parent group virial 
mass. For groups with M > 10 13 ' Mq the relative fraction of 
star forming galaxies is significantly suppressed, indicating 
that even intermediate mass environments affect star for- 
mation. Dommguez et al. (2002) presented hints toward un- 
derstanding local environment effects affecting the spectral 
types of galaxies in groups by studying the relative fractions 
of different spectral types as a function of the projected lo- 
cal galaxy density and the group-centric distance. A similar 
analysis were performed in known galaxy clusters and their 
environments in the 2dFGRS by Lewis et al (2002). 

The aim of this work is to use the Merchan & Zandi- 
varez (2002) group catalogue to obtain reliable determina- 
tions of internal and global properties of groups: luminosity 
functions of galaxies in groups, group luminosity and mass 
functions. The outline of this paper is as follows. In section 
2, we present a revised version of the 2dF Galaxy Group 
Catalogue (2dFGGC) used throughout this work. Section 3 
describes the methods and results of the luminosity func- 
tion of galaxies in groups while in section 4 description cor- 
responds to the luminosity function of galaxy groups. The 
computation of group mass function and a comparison with 
analytical models are presented in section 5. Finally, in sec- 
tion 6 we summarize our conclusions. 



2 THE 2DFGGC 

Samples of galaxies and groups used in this work are a new 
sample constructed using a revised version of the masks 
and mask software of the 2dFGRS 100k data release, which 
inclu des the /j-masks described in Colless et al. (2001) 



(see tittp://msowww. anu.edu.au/2dFGRS/Public/Release/ 
Masks/index. html). The group catalogue is obtained fol- 
lowing the same procedure as described by Merchan & 
Zandivarez (2002). In the previous work, the finder algo- 
rithm used for group identification is similar to that de- 
veloped by Huchra & Geller (1982) but modified in order 
to take into account redshift completeness and magnitude 
limit mask present on the current release of galaxies. Here 
we include the effect introduced by the magnitude complete- 
ness mask (/i-mask). In the construction of the 2dFGGC 
values of 5p/p — 80 and Vb = 200 km s _1 were used to 
maximize the group accuracy. The revised group catalogue 
comprises a total number of 2198 galaxy groups with at 
least 4 members and mean radial velocities in the range 
900 km s" 1 < V < 75000 km s" 1 . These groups have a mean 
velocity dispersion of 265 km s _ , a mean virial mass of 
9.1 x 10 13 ft" 1 M Q and a mean virial radius of 1.15 h x Mpc. 
These results show that the new identification keeps the 
mean properties of the group catalogue obtained by Merchan 
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Figure 1. bj band luminosity function of galaxies in groups from 
the 2dFGGC (arbitrary units) computed using the C~ method. 
Error bars were estimated using mock catalogues as described 
in text. Solid line shows the STY Schcchtcr fit to the data. Inset 
panel displays l<r, 2a and 3a confidence ellipses enclosing the best 
fit Schechter function parameters a = —1.13 ± 0.02 and M* — 
51og/i = -19.90 ±0.03. 



& Zandivarez (2002), as expected from the fact that the in- 
troduction of the magnitude completeness mask is a second 
order correction. 



3 LUMINOSITY FUNCTION OF GALAXIES 
IN GROUPS 

In this section we estimate the luminosity functions (LF) 
of galaxies in groups. The following analysis comprises the 
study of the whole sample of galaxies in groups and differ- 
ent subsets defined by galaxy spectral types as defined by 
Madgwick et al (2002). This classification is based on the r\ 
parameter which is very tightly correlated with the equiv- 
alent width of H a emission line, correlates well with mor- 
phology for emission line galaxies and can be interpreted as 
a measure of the relative current star-formation present in 
each galaxy. The four spectral types are defined as: 



Type 1 
Type 2 
Type 3 
Type 4 



-1.4 < Tj < 

1-1 < V < 



V < -1-4, 
1.1, 
3.5, 
Tj > 3.5. 



The Type 1 class is characterised by an old stellar popu- 
lation and strong absorption features, the Types 2 and 3 
comprise spiral galaxies with increasing star formation, fi- 
nally the Type 4 class is dominated by particularly active 
galaxies such as starburst. 



3.1 Luminosity function estimators 

In a comparative study of different LF estimators using 
Monte Carlo simulations, Willmer (1997) shows that the C~ 
method of Lynden-Bell (1971) and the STY method derived 
by Sandage, Tammann & Yahil (1979) are the best estima- 
tors to measure the shape of the LF. Moreover, Willmer's 
analysis states that the C~ is the most robust estimator, be- 
ing less affected by different values of the faint end slope of 
the Schechter parametrisation and sample size. In this work 
we use both, the C~ method to make a non parametric de- 
termination of the LF and the STY method to calculate the 
maximum likelihood Schechter fit to the LF. 

The C~ method was simplified and developed by 
Choloniewski (1987) in order to compute simultaneously the 
shape and normalisation of the luminosity function. The LF 
is obtained by differentiating the cumulative LF, ^(M). The 
function X(M) defined as the observed density of galaxies 
with absolute magnitude brighter than M, represents only 
an undersampling of the ty(M), 



dX 
X 



< 



(1) 



the key point of the method is to construct a quantity C(M), 
subsample of X(M), such that 



dX 
C 



(2) 



Linden-Bell (1971) denned that quantity as the number of 
galaxies brighter than M which could have been observed 
if their magnitude were M. Following Choloniewski (1987), 
and taking into account the sky coverage of the 2dF present 
release, the differential LF can be written as 



<<&(M)> = 



rE, 1 



Mi£[M,M+AM] 



AM 



where 



fc=2 
k 

^ = n 



C'k + w k 



Ck 



(N N 
i=l j = l 




Ci + Wi 



i+i 



(3) 



(4) 



(5) 



and R(a, S) is the redshift completeness of the parent cata- 
logue. Ck = C~ (Mfe) is denned as C(M) but excluding the 
object k itself and weighting each object by the inverse of 
the magnitude-dependent redshift completeness denned as 
w^ 1 — Cz(bj,fi) — 0.99(1 — exp(fej — fj,)) where \i = n{a,S) 
(see Norberg et al. 2002). 

In addition to the non-parametric method described 
above we also apply the STY method which uses a maximum 
likelihood technique to find the most probable parameters of 
an analytical <E>(M), in general assumed to be a Schechter, 
function: 



$(M) oc 10-°' 4(M 



-M*)(a+1) 



exp[-10-°- 4(M 



-M*), 



(6) 



The probability p, of having an object with absolute mag- 
nitude Mi is 



Pi 



*(Mi) 



j:-*(M )d M 



(7) 
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where Mb r i g ht and Mf a i nt are the brightest and faintest ab- 
solute magnitudes observable in the sample at the redshift 
of the considered galaxy with the corresponding k + e cor- 
rection. Consequently, the method seeks for the parameters 
M* and a maximizing the likelihood function 



n 



Pi- 



As we did before, we use the ^i-mask to weight each galaxy 
by l/c z (bj, jj,) to compensate for the magnitude-dependent 
incompleteness. The errors can be estimated from the error 
ellipsoid defined as 



1 2 

ln£ = ln£ max - -Xp(N), 



(9) 



where Xf)(,N) is the /3 point of the \ 2 distribution with N 
degrees of freedom. 

The comoving volumes involved in previous equations 
are estimated with the general formula 



V 



d L (z) 2 (l + z)- 2 dz 



H o J X1 y/Qo{l + z) 3 + (1 - fi - «a)(1 + zf + fl A 



(10) 



where ui is the solid angle, c is the speed of light, Ho = 
100 h km s _1 Mpc -1 and the cLl{z) is the luminosity dis- 
tance. Hereafter we adopt the cosmological model Qg = 0.3 
and SIa = 0.7. 



3.2 Luminosity function results 

In Figure [j] we show the luminosity function for all galax- 
ies in groups in arbitrary units. It should be taken into ac- 
count that the normalisation procedure is only necessary for 
the group LF. Absolute magnitudes are computed using the 
k + e mean correction for galaxies in the 2dFGRS as derived 
by Norberg et al. (2002). Error bars are estimated using 
10 mock catalogues constructed from numerical simulations 
of a cold dark matter universe according to the cosmolog- 
ical model adopted in this work with a Hubble constant 
h = 0.7 and a relative mass fluctuation as = 0.9. These 
simulations were performed using 128 3 particles in cubic co- 
moving volume of 180/i _1 Mpc per side. The solid line shows 
the STY best fit which corresponds to the Schechter param- 
eters M* = -19.90 ± 0.03 and a = -1.13 ± 0.02. Quoted 
errors are the projections of la joint error ellipse (inset plot) 
onto each axis. The inset plot of Figure |l| also shows the error 
ellipsoids defined by 2 and 3 a error dispersion. 

As pointed out by Martinez et al. (2002), the relative 
fraction of the spectral types in groups and field are different, 
being Type 1 galaxies the dominant population in groups. 
We are interested in deepening this result by studying the 
LF for the four spectral types in groups. The luminosity 
functions for different galaxy spectral types are shown in 
Figure ^. In these computations we applied a redshift cut- 
off z < 0.15 since 77 determinations are available only for 
galaxies in this redshift range. In this case, error bars were 
estimated using the bootstrap resampling technique. The 
corresponding STY best fit parameters are quoted in Table 
1, and were computed using only those galaxies brighter 
than Mb j — 5 log h = — 17 avoiding points with larger error 
bars. 

At this point, a straightforward comparison between our 



results and those obtained for field galaxies in the 2dFGRS 
by Madgwick et al. (2002) can be made. The comparison 
between the LF Schechter parameters is shown in Figure ^. 
Overall, a brightening of the characteristic magnitude is ob- 
served for galaxies in groups with a statistical significance 
~ 2a. This is particularly significant for Type 1 galaxies 
that are the main contributors to the general LF of galax- 
ies in groups (3<r). The denser environment of groups could 
be thought as the natural responsible of this brightening. In 
this particular environment, galaxy mergers seems to be the 
most probable cause since tidal interactions are more effec- 
tive when the galaxy encounter is slow. The later situation, 
is more frequently observed in groups where smaller veloc- 
ity dispersions contrast with the high velocity dispersions 
observed in rich clusters. Other processes, as ram-pressure 
(Abadi, Moore & Bower 1999) or galaxy harassment (Moore 
et al. 1996) should not be as effective as in rich cluster where 
the intra-cluster gas is much denser and the interaction rate 
is higher. 

Nevertheless, the processes mentioned above may be im- 
portant in the generation of red low mass galaxies, that 
are possibly the remnants of dynamical stripped galaxies 
in high mass systems and can be the responsible for the 
steepening of the luminosity function in clusters. The ob- 
served difference in the a parameter for Type 1 galaxies 
(lower panel of Figure ^ may be explained in this frame- 
work. In order to test this scenario, we reanalyze the LF for 
non star-forming galaxies (Type 1) splitting the sample in 
high (M > 10 14 M s ) and low (M < 10 14 M ) group mass 
subsets. A significant difference is observed at the faint end 
slope of the LF from this comparison, resulting Schechter 
a parameters of (-1.10 ± 0.06) and (-0.71 ± 0.04) for high 
and low mass subsamples respectively (Figure ^ . These re- 
sults may support our hypothesis that in groups, mergers 
are probably responsible for the M* brightening, meanwhile 
process as the ram-pressure and galaxy harassment could 
generate the increase at the LF faint end slope as observed 
for non-star forming galaxies in high mass systems. Besides 
the steepening of the faint end slope of Type 1 galaxies LF in 
higher mass systems, the galaxy LF for high and low mass 
groups do not differ appreciatively between them showing 
roughly the same behaviour as the overall LF of galaxies in 
groups (Figure [j]). 

Despite the differences between the resulting luminos- 
ity functions per spectral type, the general LF of galaxies 
in groups is quite similar to that of 2dF field galaxies ob- 
tained by Norberg et al. (2002) and Madgwick et al. (2002). 
This behaviour could be possibly due to the different relative 
abundances of Type 1 galaxies in the field and groups. In 
groups Type 1 galaxies are by far the dominant population, 
consequently, these galaxies almost determine by themselves 
M* . Meanwhile, in the field, Types 2, 3 and 4 contribute 
more strongly to the LF, generating a M* brighter than the 
corresponding to Type 1 galaxies alone. Regarding the faint 
end slope of the overall LF, the main contributors are late 
type galaxies (except for high mass groups, as stated above) 
in both, field and groups. Since there are no significant dif- 
ferences between the a parameters for late types in the field 
and groups, then the overall faint end slopes are similar. 
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Figure 2. bj band luminosity function per spectral type for galaxies in groups from the 2dFGGC (arbitrary units) computed using the 
C* _ method. Error bars are bootstrap resampling technique estimates. Solid lines show the STY Schechter fits to the data. Inset panels 
display la, 2a and 3a confidence ellipses enclosing the best fit Schechter function parameters, see Table 1. 



Table 1. STY estimates of the Schechter parameters for the luminosity functions of galaxies in groups 



Sample 


Redshift 


range 


Number of galaxies 


M* -51og(fc) 


a 


All 


0.003 - 


0.25 


14620 


-19.90 ± 0.03 


-1.13 ±0.02 


Type 1 


0.003 - 


0.15 


6202 


-19.85 ± 0.04 


-0.76 ±0.03 


Type 2 


0.003 - 


0.15 


3445 


-19.54 ±0.05 


-0.94 ±0.04 


Type 3 


0.003 - 


0.15 


1871 


-19.37 ±0.08 


-1.25 ±0.05 


Type 4 


0.003 - 


0.15 


952 


-19.40 ± 0.10 


-1.40 ±0.07 
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4 LUMINOSITY FUNCTION OF GROUPS 

The next step in this work is to analyse the luminosity distri- 
bution of groups independently of the detailed arrangement 
of luminous material within each object. In this section we 
determine individual group luminosities that allow us the 
computation of the group luminosity function and the sub- 
sequent comparison with previous determinations and pre- 
dictions from semianalytical models of galaxy formation. 



4.1 Group luminosities 

We compute group luminosities following the prescription 
by Moore, Frenk & White (1993). For each group the total 
luminosity, L t ot , is the sum of the luminosities of the galaxy 
members (Lobs) plus the integrated luminosity of galaxies 
below the magnitude limit of the survey (L C or) 



(11) 



In the computation of L COI we assume that group mem- 
bers are independently drawn from the previously computed 
Schechter fit to the luminosity function of galaxies in groups. 
Consequently, the expected luminosity of faint group mem- 
bers is 



^cor — ^Vo 



(12) 



where L Xim = W oMM e -M lim ) ^ M Um = m Mm - 25 - 
b\og{d L {z)) and M Q = 5.30 {bj band). The reliability of 
this scheme for computing luminosities of virialized systems 
was widely tested by Moore, Frenk & White (1993) using 
groups from the CfA catalogue and flux and volume limited 
mock catalogues. In order to apply this procedure to our 
group catalogue we introduce some changes that take into 
account the particular sky coverage of the parent catalogue: 
mn m is a function of the angular position of the group and 
N b B — wnere trie sum extends over the group 

members and each member is weighted with the inverse of 
the redshift completeness, Wi — 1/R(cti,5i), at its angular 
position. 



4.2 Group luminosity function 

As a first step in computing the group luminosity function 
using the C~ method, it is necessary to find out the com- 
pleteness limit in group apparent magnitude of our cata- 
logue. Rauzy (2001) proposed a new tool in order to de- 
fine this completeness limit for a magnitude-redshift sample. 
This method does not presuppose that the galaxy popula- 
tion does not evolve with time and is homogeneously dis- 
tributed in space as the V/V m ax test of Schmidt (1968). 
Rauzy's method assumes that the sample is complete in 
apparent magnitude up to a given magnitude limit mii m . 
The limiting magnitude determines a correlation between 
the absolute magnitude M and the k + e corrected distance 
modulus Z. This method is based on the definition of a ran- 
dom variable 



*(M) 



where Mn ln (Z) = mn m — Z. C, should be uniformly dis- 
tributed between and 1 and should be statistically inde- 
pendent of Z and the angular selection function. For each 
object, an unbiased estimate of ( is provided by 



= 



(14) 



TH + V 

where n is the number of objects in the sample with M < 
Mi and Z < Zi, and rii is the number of objects such that 
M < Mf im (Zi) and Z < Zi. The expectation Ei and variance 
Vi of the (i are 1/2 and (n; — l)/(12(n; + 1)) respectively. 
Then, the statistic Tc, defined as 



1/2 



(15) 



tf[M Um (Z)] 



(13) 



has mean zero and variance unity. The test consists in com- 
puting the quantity Tc on truncated subsamples according 
to an increasing apparent magnitude limit m* . While m* re- 
mains below mum the subsample is complete and Tc statis- 
tics is distributed with a mean close to zero and unity fluc- 
tuations. When m* becomes greater than m\ lra , the incom- 
pleteness introduces a lack of objects with M fainter than 
Afiim(Z), therefore Tc is expected to fall down systemati- 
cally to negative values. The behaviour of Tc for our group 
catalogue is shown in the upper panel of Figure where 
it can be seen that Tc decreases monotonously for limit- 
ing apparent magnitudes greater than mn m ~ 15.6 (vertical 
solid line), taking values under — 2(= — 2ctt c ) (horizontal 
long dashed lines) . The lower panel of the same figure shows 
the k + e corrected distance modulus Z as a function of the 
absolute magnitude for all groups, the vertical solid line cor- 
responds to the 15.6 limiting apparent magnitude cutoff as 
determined with the Tc analysis. Dashed horizontal line in 
the lower panel of Figure |^ represents a low redshift cut-off 
(z C o = 0.04) that we impose in order to avoid small volume 
undersampling effects. 

Finally, our flux-limited group sample comprises 922 
groups with the constraints 0.04 < z < 0.25 and m& , < 15.6. 
We compute for this sample the group LF with the C'~ 
method in a similar way as we did for galaxies in groups 
in Section 3, but in this case no weight is needed in the 
Cfc = C~(Mk) computation (i.e. w = 1 in equations 4 
and 5) since the choice of the magnitude cut-off ensures 
the completeness of the final sample. The resulting group 
LF is shown in Figure |(| Error bars were computed by ap- 
plying the same LF estimator to the mock catalogues in 
a similar way as explained is Section 3. In the same fig- 
ure we show the results from groups in the CfA redshift 
survey by Moore, Frenk & White (1993) and from semian- 
alytical models by Benson et al. (2000) for a A cold dark 
matter model. A general good agreement with previous re- 
sults is observed except for Mtj — 51og h > —21.5. Since the 
2dFGGC has only groups with at least 4 members, there 
is a lack of low luminosity systems that could explain the 
differences in the low luminosity tail in Figure |(S| when com- 
paring with previous works. As shown by Moore, Frenk & 
White (1993), the faint end of the group luminosity function 
is made up almost enterely of single galaxies; binaries and 
groups with three members link this to the steep bright tail 
of richer groups and clusters. We have checked whether the 
4-member criterion is affecting the faint-end of the group lu- 
minosity function by re-computing it including groups with 
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Figure 3. Comparison between the LF Schechter parameters of 
field galaxies by Madgwick ct al. (2002) (open circles) and galaxies 
in groups of this work (filled circles). The upper panel compares 
the characteristic magnitude M* — 51og(/i), meanwhile the lower 
panel is the comparison for the faint-end slope a. Both panels are 
plotted as a function of galaxy spectral type. We label as Type 
all galaxies, irrespectively of spectral type. 



Figure 5. Upper panel: the test for completeness Tq (Rauzy 
2001) applied to 2dFGGC. Apparent bj magnitudes were com- 
puted using group total luminosities following Moore, Frenk & 
White (1993) and group redshifts. A systematic decline of the 
Tq statistics can be observed beyond mij m = 15.6 (vertical line). 
Lower panel: Distance modulus Z versus absolute magnitude for 
groups. Solid line shows the apparent magnitude of completeness 
bj = 15.6. Horizontal dashed line represents the low redshift cut- 
off z > 0.04. 
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Figure 4. Luminosity functions for Type 1 galaxies in high (left 
panel) and low (right panel) mass groups samples. Solid lines are 
the best Schechter fit to the luminosity function determined with 
the STY method. 



3 and 2 members. The resulting group LF is in agreement 
with the results of Moore, Frenk & White (1993) in the 
whole absolute magnitude range. Nevertheless, it should be 
recalled that the 4 member cut-off is necessary in order to 
avoid spurious detections. As stated by Merchan & Zandi- 
varez (2002), this false detections can reach aproximately 
~ 40% when considering groups with only three members. 
Consequently, the 4-member criterion ensures a reliable de- 
termination of the luminosity function of virialised systems 
over a wide range of luminosities. 



5 MASS FUNCTION OF GROUPS 

The statistical strength of our group sample allows us an 
accurate measurement of the group mass function. In this 
section we explore the mass function of intermediate mass 
systems, using the flux-limited sample defined in the previ- 
ous section. For this purpose we use the 1/Vmax procedure, 
that is a non-parametric method consisting in weight each 
object with the inverse of the maximum comoving volume 
of the survey, V max , 



Vmax(M;; 



/■max(zi,z min ) 

I t \ dz 



dz 



(16) 



in which it remains observable given all the observational 
selection limits. The extremes of the volume integral are 
found by solving the equations 
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Figure 6. Group luminosity function (solid circles) for our flux- 
limited group sample compared with the CfA redshifts survey 
groups (open triangles) Moore, Frenk & White (1993) and the 
semianalytical models of Benson et al (2000) (dashed lines). 



Mi = 



1112 

mi 



5 log d L (z n 
5 log d L (z n 
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(17) 



that are the redshifts of objects with the same absolute mag- 
nitude, Mi, of the considered object, but with apparent mag- 
nitude at the bright and faint limits of the survey. Conse- 
quently, the differential mass function can be computed as 

dN ;(M)= £ „ ,„L (18) 



dM 



\Mi-M\<dM/2 



Vm ax (Mi)R(ai,5i) 



where Mi are the group masses and R(ai,5i) are the mean 
group redshift completeness. 

For a meaningful comparison between the resulting 
group mass function and analytical mass function predic- 
tions (Press & Schechter 1974, Sheth & Tormen 1999 and 
Jenkins et al. 2001), it is necessary to relate the group 
virial masses with the mass definition in the models. Group 
virial masses were computed inside a galaxy density contrast 
(Sp/p) g — 80 (Merchan & Zandivarez 2002), which corre- 
sponds to a mass density contrast (Sp/p) m = (5p/p) g /b, 
where b is the linear bias factor. If we estimate b = 1/<j&, as- 
suming as = 0.9, the resulting density contrast is (Sp/p) m = 
72. In order to make a comparison with the analytical predic- 
tions, the masses should be computed for a particular den- 
sity contrast depending on the cosmological model. Taken 
into account the relation between the virial density of col- 
lapsed objects in units of the critical density and the adi- 
mensional density parameter Qo, as quoted by Eke, Cole & 
Frenk (1996), we must have {8p/p) m ~ 330 for our cosmo- 
logical model, assuming spherical collapse approximation. 
To find the relation between the group virial masses and 
that corresponding to the matter overdensity 330, we as- 
sume a simple mean overdensity profile for groups that scales 



Figure 7. Group mass function for the 2dFGGC (dots). Error 
bars correspond to one cr dispersion obtained from ten mock cat- 
alogues. The solid line, dotted line and dashed line are the Press 
& Schechter (1974), Sheth & Tormen (1999) and Jenkins et al. 
(2001) prescriptions respectively. 



as r , where r is the group-centric distance. Given that 
in this case Mar then the resulting scaling relation is 
~ M 72^72/330. 
Figure [?] shows the comparison between the group mass 
function, with masses scaled using the previous relation, and 
the analytical prediction for Press & Schechter (1974), Sheth 
& Tormen (1999) and Jenkins et al. (2001). Overall, analyt- 
ical predictions are in excellent agreement with our results 
spanning two orders of magnitude in mass (10 13 /i _1 Mq <J 
M < 10 15 /i _1 M©). The differences for masses smaller than 
~ 10 13 hT 1 Mq arise as consequence of the lack of low mass 
groups, mainly due to the member number cut-off imposed 
to the 2dFGGC. 



6 CONCLUSIONS 

Here, we have applied several statistical analysis to groups 
and galaxies in groups taken from an updated version of the 
Merchan & Zandivarez (2002) group catalogue (2dFGGC). 
We have focused on an accurate determination of the lumi- 
nosity functions of galaxies in groups, and the luminosity 
and mass functions of groups taking advantage of the sta- 
tistical power of the 2dFGGC. 

In the LF computations, we have used a version of 
the Choloniewski (1987) approach to the C~ method by 
Lynden-Bell (1971), adapted to the particular sky coverage 
of the 2dFGRS. The choice of this particular estimator of 
the luminosity function is inspired in the conclusions of the 
comparative analysis of LF estimators by Willmer (1997), 
who states that the C~ method is less affected by inho- 
mogeneities in the sample. The resulting luminosity func- 
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tion for galaxies in groups (Figure |l|) is well fitted by a 
Schechter function with shape parameters M* — 5 log(ft) = 
-19.90 ± 0.03 and a = -1.13 ± 0.02 as determined by a 
STY fitting procedure. These values are quite consistent 
with those obtained by Norberg et al. (2002) for field galax- 
ies in the 2dFGRS. 

We have performed a similar analysis in subsamples of 
galaxies (Figure ^| and Table 1) defined by the spectral types 
of Madgwick et al. (2002) . In general, the characteristic mag- 
nitudes M* are shifted to higher luminosities with respect 
to the field values found by Madgwick et al (2002), irre- 
spectively of spectral type. This shift may be due to galaxy 
dynamical interactions such as mergers, which are expected 
to be much more frequent in systems with low velocity dis- 
persions as the majority of our group sample. 

The faint end slopes of the luminosity functions, a, 
are consistent with those corresponding to field galaxies ex- 
cept for low star forming, Type 1, galaxies, that show a 
steeper value in groups (Figure ^). We deepen our anal- 
ysis for Type 1 galaxies exploring the behaviour of a for 
two subsets of groups: low (M < 10 14 /i _1 Mq) and high 
(M > 10 14 ft _1 M Q ), virial masses. We observe an increase 
in the faint end slope of the Type 1 LF (a = -1.10 ± 0.06) 
for galaxies in high mass systems meanwhile for galaxies 
in low mass systems it remains closer to the global value 
(a = -0.71 ± 0.04) (Figure §). This effect could be the 
result of internal processes in higher mass system environ- 
ments such as ram-pressure and galaxy harassment, that are 
not expected to be significantly important in smaller over- 
densities. 

We have defined group luminosities following Moore, 
Frenk & White (1993), as the sum of its observed lumi- 
nosity plus a normalised integral of the galaxy luminosity 
function below the flux limit of the survey. This definition 
has proved to be reliable in the computation of group total 
luminosity using observations and artificial flux and volume- 
limited catalogues. The main aim of this computation is the 
determination of the luminosity function for our group cat- 
alogue. Since the C~ method for luminosity function esti- 
mation requires a fair selection criterion, we have used the 
Rauzy (2001) Tc statistics to determine an apparent magni- 
tude cut-off for the 2dFGGC that ensures the highest level 
of completeness. Our final flux-limited group sample com- 
prises 922 groups with apparent magnitudes brighter than 
bj — 15.6 (Figure^). The resulting group luminosity func- 
tion for this sample (Figure |^) is consistent with Moore, 
Frenk & White (1993) determination for groups in the CfA 
redshift survey and with the semianalytical model prediction 
of a A cold dark matter cosmology performed by Benson et 
al. (2000). This agreement is acceptable in the luminosity 
range Mbj < —21.5. The differences observed at fainter lu- 
minosities are mainly due to the lack of groups with less than 
4 members in the 2dFGGC. The intrinsic characteristics of 
the group finding algorithm used in the construction of the 
2dFGGC by Merchan & Zandivarez (2002), determine that 
below this limit any resulting sample has an unacceptable 
level of contamination by spurious detections. 

The flux-limited group sample adopted in the luminos- 
ity function computation is also fair enough for the determi- 
nation of the group mass function. This determination was 
achieve d by using an adapted version of the 1/Vmax that 



comparison with analytical predictions of halo abundances 
is made using a simple scheme to relate group virial masses 
and that corresponding to the adequate overdensity in our 
cosmological model. The results are displayed in Figure ^, 
where it can be seen a notorious agreement with the Press 
& Schechter (1974), Sheth & Tormen (1999) and Jenkins 
et al. (2001) mass functions for masses M > IO^^Mq. 
Again, the disagreement for low masses can be attributed to 
absence of poor groups. 

The statistical significance of the group catalogue used 
in this work allowed us to obtain very important clues about 
internal properties of intermediate mass systems and also to 
observe the level of agreement obtained from both analytical 
and semianalytical predictions in the framework of a A cold 
dark matter model. 



ACKNOWLEDGMENTS 

We are indebted with Diego G. Lambas and Hernan Muriel 
for interesting comments and discussions. We thank to Peder 
Norberg and Shaun Cole for kindly providing the software 
describing the mask of the 2dFGRS and to the 2dFGRS 
Team for having made available the current data sets of 
the sample. We would like to thank J. Coltrane. HJM, AZ 
and MEM are supported by fellowships from CONICET, 
Argentina. MJLD is supported by a fellowship from SE- 
CyT, Universidad Nacional de Cordoba, Argentina. This 
work has been partially supported by grants from Consejo 
de Investigaciones Cienti'ficas y Tecnicas de la Republica Ar- 
gentina (CONICET) , the Secretarfa de Ciencia y Tecnica de 
la Universidad Nacional de Cordoba (SeCyT), Agencia Na- 
cional de Promocion Cienti'fica de la Republica Argentina 
and Agencia Cordoba Ciencia. 



REFERENCES 

Abadi M.G., Moore B., Bower R., 1999, MNRAS, 308, 947. 
Bahcall N.A., Cen R., 1993, ApJ, 407, L49. 

Benson A.J., Cole S., Frenk C.S., Baugh CM., Lacey C.G., 2000, 

MNRAS, 311, 793. 
Biviano A., Girardi M., Giuricin G., Mardirossian F., Mezzetti 

M., 1993, ApJ, 411, L13. 
Blanton M.R., et al., 2001, AJ, 121, 2358. 
Colless M., et al. (2dFGRS Team), 2001, MNRAS, 328, 1039. 
Choloniewski J., 1987, MNRAS, 226, 273. 

Dommguez M.J., Zandivarez A. A., Martinez H.J., Merchan M.E., 

Muriel H., Lambas D.G., 2002, MNRAS, in press. 
Eke V.R., Cole S., Frenk C.S., 1996, MNRAS, 282, 263. 
Garcia A.M., 1993, A&AS, 100, 47. 

Girardi M., Borgani S., Giuricin G., Mardirossian F., Mezzetti 

M., 1998, ApJ, 506, 45. 
Girardi M., Giuricin G., 2000, ApJ, 540, 45. 

Goto T., Okamura S., Mckay T. A., Annis J., Bahcall N.A., 
Bernardi M., Brinkmann J., Gomez P.L., Hansen S., Kim 
R.S.J., Sekiguchi M.. Sheth R.K.. 2 002, to be published in 
PASJ vol.54 No.4, |astro-ph/0205413| . 

Huchra J. P., Geller M.J., 1982, ApJ, 257, 423. 

Jenkins A., Frenk C.S., White S.D.M., Colberg J.M., Cole S., 
Evrard A.E., Yoshida N., 2001, MNRAS, 321, 372. 

Lin H., Kishner R.P., Shectman S.A., Landy S.D., Oemler A., 
Tucker D.L., Schechter P.L., 1996, ApJ, 464, 60. 

Lewis I., et al. ( 2HFGRS Team), 2002, MNRASsubmitted, astro- 



conside ~s the sky coverage of the group sample. Finally, a 



ph/020333' 



© 0000 RAS, MNRAS 000, 000-000 



10 H.J. Martinez et al 



Lopez-Cruz O., Ycc H.K.C., Brown J. P., Jones C, Forman W., 

1997, ApJ, 475, 97. 
Loveday J., Peterson B.A., Efstathiou G., Maddox S.J., 1992, 

ApJ, 390, 338. 
Lyndcn-Bell D., 1971, MNRAS, 155, 95. 

Madgwick D., et al. (2dFGRS Team), 2002, MNRAS, 333, 133. 
Martinez H.J., Zandivarcz A., Dommguez M., Mcrchan M.E. & 

Lambas D.G., 2002, MNRAS, 333, 31L. 
Marzke R.O., Huchra J. P., Geller M.J., 1994, ApJ, 428, 43. 
Merchan M.E., Zandivarez A., 2002, MNRAS, in press. 
Moore B., Frenk C.S., White S.D.M., 1993, MNRAS, 261, 827. 
Moore B., Katz N., Lake G., Dressier A., Oemler A., 1996, Nature, 

379, 613. 

Marinoni C, Hudson M., Giuricin G., 2002, ApJ, 569, 91. 
Muriel, H., Valotto, C. & Lambas, D. G., 1998, ApJ, 506, 540. 
Norbcrg P., et al. (2dFGRS Team), 2001, MNRAS, 328, 64. 
Press W., Schechter P., 1974, ApJ, 187, 425. 
Ratcliffe A., Shanks T., Parker Q.A., Fong R., 1998, MNRAS, 

294, 147. 
Rauzy S., 2001, MNRAS, 324, 51. 

Sandage A., Tammann G.A., Yahil A., 1979, ApJ, 352, 82. 
Schmidt M., 1968, ApJ, 151, 393. 
Sheth R.K., Tormen G., 1999, MNRAS, 308, 119. 
Trentham N., 1997, MNRAS, 290, 334. 

Valotto C.A., Nicotra M.A., Muriel H., Lambas D.G., 1997, ApJ, 
479, 90. 

Willmer C.N.A, 1997, AJ, 114, 898. 
Zucca E. et al., 1997, A&A, 326, 477. 



© 0000 RAS, MNRAS 000, 000-000 



